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ABSTRACT 

This  paper  describes  the  application  of  a  laboratory 
based  experimental  method  [1]  to  three  dimensional 
cracked  body  problems  in  pressure  vessels  in  order  to 
determine  the  crack  shape  and  stress  intensity  fac¬ 
tor  (SIF)  distribution  along  the  crack  front  when  the 
crack  shape  is  not  known  a-priori.  Results  for  specific 
problems  are  presented  and  conditions  and  limitation 
of  the  method  are  described. 


INTRODUCTION 

Back  in  the  mid  sixties,  when  an  explosion  of  three 
dimensional  fracture  problems  solutions  developed 
along  with  the  digital  computer,  the  writer  and  his 
associates  decided  that  a  laboratory  based  experi¬ 
mental  method  should  be  developed  in  order  to  par¬ 
tially  fill  the  gap  produced  by  the  lack  of  proof  test¬ 
ing  for  many  large  structures,  such  as  ships,  bridges, 
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aircraft,  missiles  and  pressure  vessels  because  of  the 
excessive  costs  involved.  By  effecting  a  marriage  be¬ 
tween  the  equations  of  linear  elastic  fracture  mechan¬ 
ics  (LEFM)  with  frozen  stress  photoelasticity  (Ap¬ 
pendix  A),  such  a  method  has  been  developed.  Over 
the  years,  it  has  been  used  to  check  many  numeri¬ 
cal  solutions  in  a  cost  effective  manner  and,  in  the 
process  has  provided  some  novel,  or  unique  aspects. 
The  method  of  analysis  with  equations  for  all  three 
local  modes  of  stress  intensity  is  described  in  [1],  and 
is  also  presented  in  [2].  The  unique  feature  lies  in  the 
fact  that  fact  that  starter  cracks  may  be  grown  above 
critical  temperature  by  applying  service  loads  before 
cooling  the  model.  When  the  crack  reaches  its  de¬ 
sired  size  and  shape  (controlled  only  by  loads  and 
body  shape)  the  loads  are  reduced  in  order  to  stop 
crack  growth  and  stress  freezing  is  then  completed.  If 
sufficient  growth  is  permitted  to  occur,  only  Mode  I 
will  exist  along  the  crack  front.  That  is,  as  described 
by  Cotterell  [3],  [4]  the  crack  is  a  Class  I  crack,  and 
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by  using  data  sufficiently  close  to  the  crack  tip,  a  two 
term  expression  for  the  local  stresses  becomes  ade¬ 
quate  for  analysis.  The  resulting  equations  for  both 
mixed  and  pure  Mode  I  are  presented  in  Appendix 

B. 

EXPERIMENTAL  RESULTS 

Nozzle  Corner  Craeks  in  Reaetor  Vessels 

The  writer  and  his  associates  were  approached 
by  researchers  at  Delft  Technical  University  in  1976 
to  evaluate  a  test  developed  at  Delft  for  determin¬ 
ing  stress  intensity  factors  for  a  nozzle  corner  crack. 
They  used  a  fiat  plate  with  a  nozzle  located  in  its  cen¬ 
ter  with  a  starter  crack  oriented  normal  to  an  exter¬ 
nal  field  of  uniaxial  tension  [5] .  They  grew  the  crack 
in  a  steel  model  using  high  cycle  tension-tension  fa¬ 
tigue  loads.  The  result  is  shown  in  Fig  la  for  small 
and  moderately  deep  cracks  and  in  Fig  Ic  for  a  deep 
crack.  Using  monotonic  loads,  we  produced  shapes 
shown  in  Figs  lb  and  Id  which  overlaid  those  ob¬ 
tained  at  Delft. 

We  subsequently  applied  our  method  to  study 
a  similar  local  geometric  situation  where  the  nozzle 
was  located  in  a  scale  model  of  a  nuclear  pressure 
vessel.  (Fig.  2).  The  model  consisted  of  two  semi- 
cylindrical  center  sections,  each  containing  a  nozzle 
and  two  hemispherical  end  domes.  Starter  cracks 
were  oriented  as  shown  in  Fig.  3a  and  the  model  was 
loaded  by  internal  pressure.  The  crack  shapes  and 
normalized  SIF  distributions  obtained  are  shown  in 
Fig.  3b  and  Fig.  3c  respectively.  Also  shown  in  Fig. 
3b  are  the  loci  of  fitted  ellipses  used  by  analysts  in 
this  problem.  Fig.  3c  shows  how  the  flattening  of  the 
crack  shapes  near  their  mid  points  effectively  inverts 
the  SIF  distribution  along  the  crack  front.  When  a 
starter  crack  was  inserted  into  several  models  at  the 
0  =  45°  location,  (Fig.  3a)  very  different  results  oc¬ 
curred.  First,  the  planar  starter  crack  (shaded  in  Fig. 
4a)  grew  in  its  plane  near  n  to  n’.  However,  it  fol¬ 
lowed  a  gently  curved  path  from  to  to  n.  The  normal¬ 
ized  SIF  distributions  initially  showed  a  domed  dis¬ 
tribution  with  peak  values  near  ji/ j3max  ~  0.5  (Fig. 
4b).  However,  between  a' /T'  =  0.291  and  0.458,  a 
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crease  in  the  crack  surface  between  the  planar  and 
curved  parts  of  the  crack  surface  smoothed  out  and 
the  SIF  distribution  inverted  with  the  maximum  SIF 
values  occurring  at  jillimax  ~  1.0  (Fig.  4b).  It  is 
noted  here  that  the  near  tip  fringe  patterns  along 
the  curved  portion  of  the  crack  were  not  symmetri¬ 
cal  with  respect  to  the  crack  surface,  but  no  rotation 
of  the  fringe  loops  was  present  to  indicate  a  shear 
mode,  the  lack  of  symmetry  in  the  fringe  loops  was 
due  to  a  varying  non-singular  tensile  stress  parallel 
to  the  crack  growth  direction  and  was  not  indicative 
of  a  shear  mode  [6].  Generally,  the  shape  taken  by 
the  growing  crack  is  the  one  which  minimizes  the 
gradient  in  the  SIF  along  the  flaw  border. 

Craeks  in  Roeket  Grain 

In  all  of  the  cases  discussed  above  the  starter 
cracks  were,  in  fact.  Class  I  Cracks.  Since,  how¬ 
ever,  they  are  introduced  by  striking  a  sharp  blade 
held  normal  to  a  surface,  and  the  resulting  crack  em¬ 
anates  from  the  blade  tip,  the  loading  for  the  starter 
crack  is  a  dynamic  wedging  force  and  the  induced 
stress  field  may  differ  from  internal  pressure.  This 
might  occur  if  one  drives  the  starter  crack  too  far  in. 
If  the  load  is  then  changed  to  internal  pressure,  the 
new  load  may  produce  a  shear  mode  at  the  starter 
crack  tip.  In  such  cases,  the  crack  will  then  turn 
in  order  to  eliminate  the  shear  mode  as  it  grows. 
An  example  of  this  type  of  behavior  follows.  Figure 
5  shows  a  1/10  scale  model  of  the  crossection  of  a 
cylindrical  motor  grain  model.  Two  off  axis  (Fig. 
5)  cracks  were  introduced  in  different  models  but  in 
identical  loci  of  each  of  the  models.  Plan  and  cen¬ 
ter  section  SS  views  of  the  cracks  are  shown  in  Fig. 
6.  The  crack  in  Model  4  has  not  grown  as  far  as 
the  one  in  Model  8.  If  one  assumes  that  the  crack 
in  Model  8  represents  a  later  stage  in  the  growth  of 
the  crack  in  Model  4,  one  can  construct  the  follow¬ 
ing  scenario.  Table  I  gives  pertinent  data  from  the 
frozen  stress  tests  on  the  two  cracks.  The  normalized 
(SIFs)  Fi  were  computed  using  an  approximation  of 
an  elliptical  integral  of  the  second  kind  ^/Q,  which 
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amounts  to  assuming  that  the  cracks  were  planar 
semi-elliptic  shapes.  While  this  assumption  does  in 
fact  hold  for  symmetric  cracks  (Fig.  5)  which  en¬ 
joy  both  geometric  and  load  symmetry,  such  is  not 
the  case  of  the  “off  axis”  cracks  pictured  in  Fig.  6, 
as  noted  by  the  direction  change  in  the  crack  fronts 
shown  in  the  section  SS  cuts.  The  influence  of  these 
direction  changes,  while  accommodated  in  the  ex¬ 
perimental  data,  was  ignored  in  computing  the  val¬ 
ues  of  for  use  in  determining  F,.  The  form  of 
^/Q  was  developed  by  Newman  [7].  We  note  that, 
in  our  scenario.  Model  4  represents  an  earlier  state 
in  the  growth  of  the  off-axis  crack  and  that  Model 
8  indicates  a  larger  crack.  Model  4  exhibited  both 
Fi  and  F2  values  all  along  the  crack  front  except  at 
the  fin  boundary  where  only  Fi  existed,  and  section 
SS  shows  the  crack  front  turning.  In  Model  8,  how¬ 
ever,  we  see  a  series  of  river  marks  on  both  sides  of 
the  crack  suggesting  the  presence  of  Mode  III.  No 
data  were  taken  in  these  regions  however,  since  it 
was  conjectured  that  the  limitations  of  linear  elastic 
fracture  mechanics  may  have  been  exceeded.  At  all 
other  locations  along  the  Model  8  crack  front,  how¬ 
ever,  only  pure  Mode  I  was  measured  suggesting  that 
the  crack  had  grown  sufficiently  to  eliminate  all  shear 
modes.  The  section  SS  in  Model  8  suggests  that  the 
turn  in  Model  4  (which  was  accompanied  by  Mode 
II)  has  been  completed.  This  means  that  the  crack 
in  Model  8  has  achieved,  or  nearly  achieved.  Class 
I  status  with  turning  completed  and  will  now  pro¬ 
ceed  in  the  new  direction  indicated  by  Section  SS 
(Fig.  6).  It  should  also  be  noted  that  in  Model  8,  a 
slight  dimple  appears  in  the  crack  front  in  the  areas 
of  the  river  markings,  suggesting  that  the  effect  of 
shear  modes  is  to  retard  temporarily  the  growth  of 
the  crack. 

Such  a  scenario  as  described  above  leads  to 
an  interesting  conclusion.  Both  frozen  stress  exper¬ 
iments  and  finite  element  analysis  point  to  the  fin 
corner  as  the  highest  stressed  point  in  the  uncracked 
loaded  model.  However,  any  crack  initiating  there 
will  contain  shear  modes.  On  the  other  hand,  a  crack 
initiating  on  the  fin  axis  (i.e.  the  symmetric  crack  in 
Fig.  5)  will  show  no  evidence  of  shear  modes  (Fig. 
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7)  and  is  a  Class  1  crack  when  it  appears.  As  such, 
one  expects  such  a  crack  to  be  far  more  mobile  and 
dangerous  than  one  emanating  from  the  fin  corner. 
Thus,  one  cannot  predict  a  fracture  on  the  basis  of 
uncracked  model  stress  analysis  alone.  According 
to  Cotterell,  all  Class  II  cracks  (ie  here  those  un¬ 
der  shear  modes)  eventually  become  Class  I  cracks 
and  this  concept  is  borne  out  by  these  results.  The 
symmetric  crack  shown  in  Fig.  5  is  already  a  Class 
I  crack  due  to  absence  of  shear  modes  and  will  not 
turn  during  growth. 

On  the  basis  of  the  examples  presented  above 
(and  other  cases  studied),  we  conclude  that  stable 
crack  growth  in  pressure  vessels  proceeds  under  pure 
Mode  I.  Moreover,  cracks  grown  above  critical  tem¬ 
perature  during  the  stress  freezing  process  in  pho¬ 
toelastic  models  appear  to  overlay  those  produced  in 
steel  structures  by  high  cycle  tension-tension  fatigue 
when  test  scales  are  identical. 


Table  1 


Loads' 

Crack  Description^ 
(dimensions  in  mm) 

T, 

II 

Vna  i^  1,2 

Depth  (d) 

Surface  (s) 

Model  4 

Top 

Bottom 

Fin  tip  surface 

Off-axis  inclined 

P  =  88.97  N 

P„„,  =  0.049  MPa 
Psf=  0.035  MPa 

a  =  8.71  Aa  =  2.18 

c=  11.15  Ac  =  3.02 

a/c  =  0.78  a/t  =  0.23 

F|  =  2.09 

Fi  =  0.53 

^s20° 

2.88 

ns  20° 
2.81 

/  (wp> 

/i* 

//l.OA 

lid  \/  \ 

P  =  88.97  N 

Po„,  =  0.013  MPa 

Models 

Off-axis  inclined 
a=12.50  Aa  =  3.4 

c  =  2l.l  Ac=10.4 

2.19 

ns  30° 

2.33 

ns  30° 

2.34 

SsaxU/ 

T* 

(botion) 

Psf  =0.049  MPa 

a/c  =  0.59  a/t  =  0.34 

*  Fi  values 

Notations: 


1 .  P  =  axial  compressive  load 

Pmax  =  maximum  internal  pressure  to  grow  crack 
Psf  =  stress  freezing  pressure 

2.  a  =  crack  depth;  Aa  =  crack  length  growth 

c  =  half  length  of  crack  in  fin  tip  surface;  Ac  =  half  crack  growth  in  fin  tip  surface 

3.  -v/Q  =  approximation  of  elliptic  integral  of  second  kind 


/  \1.65 

(  a  ] 

a 

/  nI.s; 
/  CL  \ 

a 

-<1  0=1.464  - 

->1 

UJ 

c 

\CJ 

c 

All  fiaws  were  characterized  as  semi-elliptic  flaws  of  depth  a  and  length  2c. 

However,  off-axis  cracks  were  neither  perfectly  semi-elliptic  nor  planar. 

4.  ns  =  near  surface  slices  were  used  in  place  of  surface  slices  to  avoid  a  possible  surface 
fringe 

5.  Slices  at  0  =  60°  revealed  a  shear  mode  in  Model  4  but  not  in  Model  8. 
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APPENDIX  A-  Frozen  Stress 
Photoelasticity 

When  a  transparent  model  is  placed  in  a  circu¬ 
larly  polarized  monochromatic  light  field  and  loaded, 
dark  fringes  will  appear  which  are  proportional  to  the 
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applied  load.  These  fringes  are  called  stress  fringes 
or  isochromatics,  and  the  magnitude  of  the  maxi¬ 
mum  in-plane  shear  stress  is  a  constant  along  a  given 
fringe. 

Some  transparent  materials  exhibit  mechanical 
diphase  characteristics  above  a  certain  temperature, 
called  the  critical  temperature  (T^).  The  material, 
while  still  perfectly  elastic  will  exhibit  a  fringe  sen¬ 
sitivity  of  about  twenty  times  the  value  obtained  at 
room  temperature,  and  its  modulus  of  elasticity  will 
be  reduced  to  about  one  six-hundredth  of  its  room 
temperature  value.  By  raising  the  model  tempera¬ 
ture  above  T^,  loading,  and  then  cooling  slowly  to 
room  temperature,  the  stress  fringes  associated  with 
Tc  will  be  retained  when  the  material  is  returned  to 
room  temperature.  Since  the  material  is  so  much 
more  sensitive  to  fringe  generation  above  than 
at  room  temperature,  fringe  recovery  at  room  tem¬ 
perature  upon  unloading  is  negligible.  The  model 
may  then  be  sliced  without  disturbing  the  “frozen 
in”  fringe  pattern  and  analyzed  as  a  two-dimensional 
model  but  containing  the  three-dimensional  effects. 
In  the  use  of  the  method  to  make  measurements  near 
crack  tips,  due  to  the  need  to  reduce  loads  above 
critical  temperature  to  preclude  large  local  deforma¬ 
tions,  and  the  use  of  thin  slices,  few  stress  fringes  are 
available  by  standard  procedures.  To  overcome  this 
obstacle,  a  refined  polariscope  is  employed  to  allow 
the  tandem  use  of  the  Post  and  Tardy  methods  to 
increase  the  number  of  fringes  available  locally. 

In  fringe  photographs,  integral  fringes  are  dark  in 
a  dark  field  and  bright  in  a  bright  field. 

Appendix  B 

Mode  I  Algorithm 

Beginning  with  the  Griffith-Irwin  Equations,  we 
may  write,  for  Mode  I,  for  the  homogeneous  case, 

fa (^)  {i,j  =n,z) 

{ZTTT)  2 

(1) 

where: 

(Ty  are  components  of  stress 
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Ki  is  SIF 

r,9  are  measured  from  crack  tip  (Fig.  B-1) 
(T°.  are  non-singular  stress  components 
Then,  along  6  =  tt  12,  after  truncating  (T,j 


_  Ffi  ,  Kap 

—  - ; — r  T  —  - - 


(2) 


where: 

r°  =  and  is  constant  over  the  data  range 

Kap  =  apparent  SIF 

=  maximum  shear  stress  in  nz  plane 


Normalizing  with  respect  to  a, 


.!l™o|  ^  {Ki,K2,rm,9m  ,Ti  j  )  |  =0 

(4) 

which  leads  to: 


4 

3 


cot  20^ 


(5) 


By  measuring  0^  which  is  approximately  in  the  di¬ 
rection  of  the  applied  load,  K2IK1  can  be  deter¬ 
mined. 

Then  writing  the  stress  optic  law  as: 


=  ^  +  (3) 

(T(7raj2  (T(7raj2  o  \a  j 

where  (Fig.  B-1)  a  =  crack  length,  and  a  =  remote 

normal  stress 

Kap  pr  . 

i.e  - -  vs.  W-  IS  linear. 

<7(710)2  V  O' 


From  the  Stress-Optic  Law,  =  n//2t  where, 
n  =  stress  fringe  order, 

/  =  material  fringe  value,  and 
t  =  specimen  (or  slice)  thickness 
then  from  Eq.  2 


max  ^ AP 

2t  (Sirr)  ^ 

where  K*^  is  the  mixed  mode  SIF,  one  may  plot 
vs.  ^/r|a  as  before,  locate  a  linear  zone  and 

(j(7ra)  2 

extrapolate  to  r  =  0  to  obtain  K* .  Knowing,  K* , 
K2IK1  and  6'^,  values  of  Ki  and  K2  may  be  deter¬ 
mined  since 


K*  =[{Kisin0l,+2K2Cos0l,f 

+  {K2sin0pf]^  (6) 


IL.4P=Cr(87rr)^=^(87rr)^ 

where  Kap  (through  a  measure  of  n)  and  r  become 
the  measured  quantities  from  the  stress  fringe  pat¬ 
tern  at  different  points  in  the  pattern. 

In  part  of  the  present  study,  instead  of  normal¬ 
izing  K  with  respect  to  <f(7ra)^/^,  we  have  selected 
Psf\/'^a/Q  as  the  normalizing  factor  where  is 
an  elliptic  integral  of  the  second  kind  approximated 
here,  as  shown  in  Table  I.  An  example  of  the  deter¬ 
mination  of  Fi  in  Table  I  from  test  data  is  given  in 
Fig.  B-2. 

Mixed  Mode  Algorithm 

The  mixed  mode  algorithm  was  developed  (see 
Fig.  B-3)  by  requiring  that: 


Knowing  K*  and  6!^,  Ki  &  K2  can  be  determined 
from  Eqs.  5  and  6.  Details  are  found  in  Smith  and 
Kobayashi  [2]. 


Fig.  B-1  Near  Tip  Notation  for  Mode  I. 


Approved  for  public  release;  distribution  unlimited. 


5 


Fig.  1  Comparison  of  Crack  Shapes  from  Delft  Fa¬ 
tigued  Steel  Model  with  those  grown  mono- 
tonically  in  photoelastic  model. 


Fig.  2  Photoelastic  Model  of  Nuclear  Pressure  Ves¬ 
sel  Containing  Two  Diametrically  Opposite 
Nozzles. 
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VESSEL  WALL 


0^*  6.74  mm  ( Test  H-A  )  (c) 


a) 


nv  SUCE 

normal  to  crack  border 


Fig.  3  a)  Starter  Crack  Locations  in  Nozzle  for  Thin 
Walled  Pressure  Vessel  b)  Flattening  of  Crack 
Shapes  with  Increasing  Depth  for  6  =  0°  c) 
Corresponding  SIF  Distributions  p  is  internal 
pressure) . 


Fig.  4  a)  Artist’s  sketch  of  planar  crack  initiated  at 
6  =  45°  b)  SIF  Distributions  for  increasing 
crack  depths. 
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length  of  cylinder  376  mm 


Fig.  5  Cross-section  of  scaled  photoelastic  model  of 
motor  grain  showing  critical  crack  locations. 


Section  S  -  S 


Starter  crack  F$ 


Mixed  mode  region 

magniflcatinn  factor:  1.25 
Model  4  Off-Axis  Inclined  Crack  Showing  Starter  Crack  and  Final  Crack  Front 


Eliminating  Mode  III 

magniiication  factor:  1.25 

Model  8  Off-Axis  Inclinetl  Crack  Showing  Starter  Crack  and  Final  Mode  I  Crack  Front 


Fig. 


FS  -  fin  surface 
<'  -  crack  front 

D  -  camera  view  of  the  photograph 


6  Plan  view  of  off-axis  crack  profiles  and  path 
of  crack  at  it’s  center  (SS). 


MDZ  -  IVhterial  Damage  Zone 
EVC  -  Edge Vie wofCmck 

Fig.  7  Plan  and  crossectional  view  of  a  symmetric 
crack. 
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